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1. Motivation, background, and objectives 

After publication of the Kolmogorov (1962) refined similarity hypotheses, the 
small-scale intermittency of the energy dissipation field became a central problem 
in fully developed turbulence (FDT). This phenomena has been studied in many 
different ways, e.g. by searching for corrections to scaling exponents in the inertial 
range velocity structure functions (see reviews in Monin & Yaglom 1975, Kuznetsov 
& Sabelnikov 1990). A direct measure of this intermittency is, however, available 
by studying the local rate of energy dissipation, and it may be quantitatively char- 
acterized by the intermittency exponent /x (Nelkin 1981). 

The first description of the intermittent field was proposed by Kolmogorov (1962), 
who introduced the energy dissipation e r averaged over a segment r: 


r 

£r( x ) = ~ + C 1 ) 

0 

where x is the longitudinal coordinate, and the local value of the energy dissipation 
e(x) is defined in the standard way as 




du j \ 

dx,J ’ 


( 2 ) 


where v is the kinematic viscosity, and tt,, i = 1, 2, 3, denote the velocity compo- 
nents in the directions x t . To analyze experimental data and compare them with 
numerical modeling, the one-dimensional sections of the energy dissipation field in 
the x = xi direction will be considered throughout this paper. 

It was assumed by Kolmogorov that e r has the normal probability density dis- 
tribution (pdd) if r belongs to the inertial range, i.e., r\ <C r <C L where L and rj 
denote the integral and Kolmogorov viscous scales, respectively. This assumption 
is known as the log-normal model. Kolmogorov further assumed that the variation 
of ln£ r obeys a logarithmic scaling 


crf ntr - A + Hh In (L/r), rj < r < £, (3) 

where A is some function of the local flow conditions, and /i* is the intermittency 
exponent which was assumed to be a universal constant (if there is no intermittency, 
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fik = 0). Different subscripts to fi are adopted to identify different definitions of 
the intermit tency exponent. 

Novikov &: Stewart (1964) showed that 

< > oc r”'* 4 * , 7] <C r <C L. (4) 

Hereafter the angular brackets denote average over x. Novikov (1990) demonstrated 
that under some weak assumptions /j* = fi € . 

It was also found (Monin &: Yaglom 1975) that the correlation function of energy 
dissipation i? £ (r) obeys the power-law scaling at separations r within the inertial 
range 


R e {r) = < e(x) e(x + r) > = C c < e > 2 ( r/L ) /ir , r] < r < I, (5) 

where C e is assumed to be a universal constant, and fi r = fi € . 

Monin Yaglom (1975) further assumed that at very high Reynolds numbers, 
where fluctuations of e(x) are much larger than < e >, 


R c {r) « 5 e (r) = < [e(x)~ < e >} [e(x + r)- <£>]>oc r tj < r < X, 

(6) 

and the energy spectrum of e (which is the Fourier transform of B e ) should behave 
as 

oo 

E e (k) — — J B € (r) cos (kr) dr cx 1/L k l/rj, (7) 

o 

where k is the wave-number. 

Relations (3)-(7) are commonly used to estimate fi. Several more methods for 
such an estimate are also well-known, e.g. by using the six-order velocity structure 
function (Monin &; Yaglom 1975), the breakdown coefficients (Novikov 1990), etc. 
but in this brief paper only methods based on (3)-(7) are considered. 

A diverse body of measurements of /i in different laboratory flows as well as 
in the atmosphere and ocean has been reported during the last three decades (for 
comprehensive review see Kuznetsov & Sabelnikov 1990, Gibson 1991, Nelkin 1994). 
The reported values of ft vary from 0.15 to 0.7, and this scatter certainly exceeds 
the measurement errors. Some theoretical considerations as well as the scatter in 
experimental results caused a doubt about universality of the exponent /i (e.g., 
Kraichnan 1974, Nelkin 1994). However, as far as we know, nobody has posed an 
obvious question: Is the intermit tency exponent /i a unique constant, i.e., are the 
values ji/t, fi r , ^6, and fi t the same at high Reynolds numbers, or do they create 

a set of different (and perhaps independent) exponents? This paper addresses the 
above question using the high Reynolds number experiments. 

The second objective of this paper becomes clear from the following considera- 
tions. It is commonly assumed (e.g., Monin & Yaglom 1975, Nelkin 1981) that the 
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longitudinal component du/dx should contain most of the important dynamical in- 
formation on the dissipation field. In other words, it is assumed that the exponents 
p defined by Eqs. (3)-(7) would not be changed if the true local rate of energy 
dissipation e(x) were replaced by the one-dimensional surrogate 

e(x) = 15 v {du/dx) 2 . (8) 

Nelkin (1981) commented, “this rather weak assumption seems plausible, but is 
difficult to test experimentally”. We are not so optimistic about the weakness 
of this assumption. Indeed, it is well-known (e.g., Hosokawa and Yamamoto 1991) 
that pdd of e(x) and £(x) are qualitatively different. The second statement, namely 
“difficult to test experimentally”, is completely valid. In most experiments nothing 
more than u(x) has been recorded, and in some experiments simultaneous measure- 
ments of u(x) and u 2 (x) were executed. There are also a few experiments at low 
Reynolds numbers (see Tsinober et al. 1992, and references therein) where all terms 
in Eq. (2) were directly measured with multi-wire probes. However, the reliability 
of such measurements, especially at high Reynolds numbers, still seems to be rather 

questionable. . 

The remarkable success of direct numerical simulations (DNS) of FDT during the 
last decades offers another approach to such a test. In spite of well-known limitations 
(relatively low Reynolds numbers, somewhat arbitrary boundary conditions, etc.), 
DNS provides an exact solution of the Navier-Stokes equations and gives complete 
information on the three-dimensional flow field. Using the numerical data base, 
one can estimate (at least qualitatively) how adequately the true energy dissipation 
e(x) is represented by the one-dimensional surrogate e(x). To execute such a test is 
a second objective of this paper. This approach contradicts traditional belief that 
measurements give the “final truth”, and any theoretical or numerical result should 
be tested by comparison with experiment. We believe that numerical modeling (now 
fairly called numerical experiment) becomes so powerful and reliable that it can be 
used to test physical experiments. First results from this approach are reported in 
this paper. 

2. Accomplishments 

2.1 Experimental and numerical data bases 
The experimental data base includes measurements in the atmospheric surface 
layer and in a large wind tunnel at very high Reynolds numbers. 

The experiment in the atmospheric surface layer was executed by Dr. Steven 
Oncley from the National Center for Atmospheric Research (Oncley 1992). After 
preliminary processing, six time series were chosen from the total record (Praskovsky 
& Oncley 1994a). Four of these series are analyzed in the present work. The second 
experiment was executed in the large wind tunnel of the Central Aerohydrodynamic 
Institute (Moscow, Russia). Measurements have been done in the mixing layer and 
in the return channel of the wind tunnel (Praskovsky et al. 1993). 

The main flow characteristics of the experimental data base are listed in Table 1 
where the abbreviations ML, RC, and ASL denote the mixing layer, return channel, 
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and atmospheric surface layer, respectively, and numerals after ASL correspond to 
the sequence of the time series (in accordance with that in Praskovsky & Oncley 
1994a). U is the mean longitudinal velocity, and a $ denotes the rms value of any 
quantity <p. The Taylor A and Kolmogorov ;/ scales, and the Reynolds number R\. 
are defined with standard formulas: A = 0 u l°du/di, = (v 3 / < £ >) 1/4 , and 
R\ = o u \jv. Other quantities in Table 1 will be defined later. Taylor’s hypothesis 
was used to convert from temporal to spatial coordinates. 


Time series 

ML 

RC 

ASL-2 

ASL-3 

ASL-5 

ASL-6 

U , m/s 

7.87 

10.8 

6.58 

8.10 

12.9 

14.5 

<7 U , m/s 

1.67 

1.03 

0.693 

1.10 

1.82 

2.08 

L , m 

1.3 

4.8 

42 

51 

99 

77 

< £ >, m 2 /s 3 

1.90 

0.115 

0.0235 

0.0322 

0.140 

0.128 

A, cm 

1.8 

4.6 

6.5 

9.0 

7.0 

8.3 

R\ x 10- 3 

2.0 

3.2 

3.3 

6.9 

9.2 

12.7 

t ], mm 

0.21 

0.41 

0.58 

0.55 

0.37 

0.37 

P7 ± <*** 

0.20±0.01 

0.19±0.01 

0.20±0.01 

0.20±0.01 

0.20±0.02 

0.20±0.01 

Wb ± 

0.38±0.01 

0.32±0.01 

0.35±0.02 

0.29±0.03 

0.33±0.02 

0.30±0.02 


0.59±0.03 

0.60±0.02 

0.57±0.04 

0.61±0.02 

0.56±0.02 

0.55±0.01 

]±k ±<?n k 

0.40±0.01 

0.22±0.03 

0.19±0.04 

0.22±0.02 

0.22±0.02 

0.23±0.02 

i3 ±crj 

0.26±0.01 

0.24±0.01 

0.23±0.02 

0.19±0.01 

0.18±0.01 

0.14±0.01 

77 ± 

0.42±0.01 

0.38±0.01 

0.41±0.02 

0.38±0.01 

0.36±0.01 

0.35±0.01 


Table 1. Main turbulence characteristics of analyzed time series. 

The DNS of homogeneous, isotropic, equilibrium flow fields were executed by 
Dr. Robert Rogallo from the NASA Ames Research Center. A description of this 
data base, covering a range of R\ from 35 to 168 can be found in Jimenez et al. 
(1993). Two simulations at R\ = 94 and 168 are analyzed in the present work. 

2.2 Experimental results for £(x) 

2.2.1 Correlation functions and spectra 

Two examples of the correlation functions R e (r) and B c (r) at the lowest and high- 
est R\ are presented in Fig. 1 . The local values of scaling exponents are estimated 
using the logarithmic derivatives: 


<*[lgfl e (r)] 

d[lgr] 


Vb{r) = 


d[lgr] 


He(k) = 1 + 


<*[lg £«(*)] 
<%*] ' 


(9) 


It follows from Eqs. (5)-(7) that p r , fx b , and fi t should be constants within the 
inertial range. The measured values of p r (r) and pfc(r) are presented in Fig. 2. 
The constant-value regions of these exponents are ill-defined, especially for n b (r). 
However, for each time series one can choose a range of separations r where the 
exponents are approximately constant. Within these limits, deviations of pt from 
constant can be attributed to measurement uncertainty. The mean values JR and 
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FIGURE 1. Correlation functions of the energy dissipation. •, R € (r)\ o, B e (r). 



r lr\ r/r| 

FIGURE 2. The local values of scaling exponents ^ r (r) and fib( r )- Note the 
expanded linear scale of the ordinate at this and other similar figures. * , ML; 
• , RC; o , ASL-2; □ , ASL-3; a , ASL-5; o , ASL-6. 


standard deviations cr^., (i = r, b) were estimated by averaging over these ranges, 
and the results are listed in Table 1. It is seen that /Jr is significantly smaller than 
Jib . This result is not surprising. Anselmet et ai (1984) found fi r = 0.18 and fib = 
0.48 in the round jet at R\ = 536. Kuznetsov et al. (1992) found fi r » 0.14, fib ~ 
0.45, and /i e « 0.60 in the mixing layer at R\ « 1700. The equality R e (r) % B e (r) 
is based on the assumption that both J? e (r) and B e (r ) are much larger than < e > 2 
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k n 


FIGURE 3. Spectra of the energy dissipation, and the local values of scaling 
exponent fi e (k). For symbols see Fig. 2. 


in the inertial range (Monin &: Yaglom 1975). This assumption is not valid at any 
finite Reynolds number. It is seen in Fig. 1 that at big enough r, say at r/rj & 10 3 , 
B € (r) is smaller than < e > 2 . Hence one should always obtain pt r < There is 
a contradiction, recently noted by Praskovsky & Oncley (1994b). The correlation 
functions are connected by the exact relation: 

B e (r) = R e (r) - < £ > 2 (10) 

It obviously follows from this relation that at any finite R\ where < e > 2 is of 
the order of (or even larger than) B e (r) at rj r T, only one of the functions, 
R e (r) or 5 c (r), may obey the power-law scaling. As was mentioned above, the 
constant-value regions for both ^ r (r) and /i&(r) are ill-defined. Taking into consid- 
eration Eq. (10), one can suggest that the present experiments as well as all previous 
measurements do not prove power-law behavior of any correlation function. The 
results demonstrate that the power laws (5) and (6) reasonably approximate both 
correlators within the inertial range. 

Typical energy spectra E € (k ) and local values of the scaling exponents fi e (k) 
are presented in Fig. 3. It is seen that at large enough k the power law (7) gives 
reasonable approximation of E e (k). The mean values /17 and the standard deviation 
cr Me are estimated over the ranges where fi e (k) are approximately constant, and the 
results are listed in Table 1. It is seen that J3S7 are roughly twice 7*6- This result 
agrees with previous studies, and it is explained by Nelkin (1981). It follows from 
Eq. (7) that if scaling (6) is valid for all r from 0 to infinity. In reality, 

B c ( 0)/ < e > 2 increases with no limit as R\ — * oo while the inertial range values 
of B c (r)/ < £ > 2 increase more slowly. In other words, B e (r) has a strong peak 
at r < 7 ) y and this peak increases as R\ increases. As was shown by Kuznetsov 





275 


Experimental & numerical study of the intermittency exponent fi 



FIGURE 4. Scaling exponents averaged over the inertial range. 

Error bars at this and other similar figures correspond to the rms value of the 
exponent. o , /Jr* a , Jib ; o , /I7- 

et al (1990), fib is always smaller than /i e , and this difference may increase when 
R\ increases. 

The values of JTJ, 7*6 > and 7*7 are plotted in Fig. 4. Beyond any doubt, /i r , 
and /i e are different exponents at R\ up to 12,700. Variation of these exponents 
with R\ and the type of flow is quite small. As a first approximation, one can 
consider the exponents fi r m 0.2, fib « 0.3, and /i e « 0.55 to be universal (but 
different!) constants at the highest Reynolds numbers currently attainable, say at 
R\ > 2xl0 3 . At Reynolds numbers too high to be attained on this planet, the 
different /i’s could still be the same. 

2.2.2 Test of the log-normal model 

The log-normal model includes the statement that ln£ r has a normal distribution 
for r/ < r < I (Kolmogorov 1962). Note that this assumption and Eq. (3) are two 
independent hypotheses, i.e., the log-normal model may be valid while Eq. (3) is 
not, and vice-versa (see Castaing et al. 1990 for more explanation). 

The log-normality of lne r was tested in numerous experiments (see reviews in 
Monin & Yaglom 1975, Gibson 1991), and this was found to give a reasonable 
approximation of reality in different turbulent flows. The present experiments give 
similar results. As an illustration, several distributions for segments r from the 
inertial range are presented in Fig. 5. It is surprising that deviations of experimental 
results from the log-normal model increase as R\ increases. 

Eq. (3) was also reported by many authors to be valid (see Gibson 1991 and 
references therein). However, Castaing et al. (1990) presented experimental and 
theoretical evidences that at small r (up to the viscous range) the variation of In s r 
obeys the power- law scaling 
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(lne r - < In £ f > ) / o |nt . ( In e f - < In e r > ) / ^ 


FIGURE 5. Pdd of 111 e T . Solid lines correspond to the normal distribution. 
RC: o , r/i) = 23; A , 130; o , 927; ASL-6: a,r/r) = 33; a , 180; o , 1025. 


2 

a \m 


cx r 




(ii) 


Hypotheses (3) and (11) are tested in Fig. 6. The local values of /xjt(r) and /3(r) 
are defined as 


Hk(r) 


±feii s,r) = - 

d[ln r] ’ d[lg r] 


( 12 ) 


It is seen that at large r, say r/r; > 200, the Kolmogorov (1962) assumption (3) 
provides a good approximation of the experimental results. At small r Eq. (3) does 
not agree with the present experiments while the power-law scaling (11) works quite 
well. The values of JT^ and (3 averaged over regions where they are approximately 
constant are listed in Table 1 and presented in Fig. 7. One can see that /x* does 
not depend on R\ and the type of flow at R\ > 3 x 10 3 , and is approximately 
equal to 0.22. The exponent j3 in (11) does not reveal any dependence on the type 
of flow (measurements by Castaing et ai 1990 in the round jet at R\ = 852 and 
in the Modane wind tunnel at R\ = 2720 are also presented in Fig. 7). However, 
this exponent strongly depends on R\, and the dependence may be approximated 
as j3 oc which is different from the relation /? oc 1/ln R\ suggested by 

Castaing et ai (1990). 


2.2. 3 Measurements of the second-order moments 

The measured values of < £ri x ) > 35 we ^ as the local values of scaling exponent 
fi e (r) — — d[ lg < € 2 r >]/d[lnr] are presented in Fig. 8. No region of power-law scaling 
is seen within the inertial range; i.e., /x ff (r) decreases monotonically over the range 
of r/r) from say 5 to 20,000. Contradiction of this result to all previous studies is 
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FIGURE 6. Variation of lne r , and the local values of /i*(r) and /?(r). 


apparent. No previous papers reported the local values /i c (r), but only log- log plots 
such as that in Fig. 8. Such plots may be deceptive; a range of slow change in fi £ (r) 
may be erroneously considered as implying fi e (r) « constant. 

It was recently suggested by LVov and Procaccia (1994) that scaling (4) should 
be applied to the centered second-order moment, i.e., the relation 

< [£ r (x)— < e >] > 2 <x r“ 7 ‘ , 7] <C r <C L (13) 

should be used instead of (4). This suggestion is tested in Fig. 9. It is seen that 
7 ff (r) = d[lg< (e r — < £ >) 2 >]/d[lnr] does reveal some range of being approxi- 
mately constant. The values of 77 averaged over these ranges are presented in 
Fig. 10. It is seen that 7* slightly depends on Reynolds number. 
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FIGURE 7. Averaged values of the exponents fi * and j3. corresponds to 

l/3 

f3 oc R x 7 ; o , the present results; A , Castaing et ai (1990). 



FIGURE 8. The second-order moments < €^(x) > and the local values of scaling 
exponent For symbols see Fig. 2. 


2.3 Analysis of the DNS data base 

Two DNS of isotropic homogeneous turbulence at R\ = 94 and 168 (Jimenez et 
al. 1993) are analyzed. The purpose of this analysis is a comparison of statistical 
characteristics of the “true” energy dissipation field £(x), Eq. (2), and its one- 
dimensional surrogate e(x), Eq. (8). Significance of such a comparison is illustrated 
in Fig. 11 where pdd P(e) of e(x) and P(i) of i(x) are presented. One can see a 
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r/r) r / rj 


FIGURE 9. The centered second-order moments < [t r (x)— < e r >] 2 >, and the 
local values of scaling exponent 7 «(r). For symbols see Fig. 2. 



FIGURE 10. Exponent 7 * averaged over the inertial range. 

qualitative difference of these functions: P(e) — ♦ 0 as e — ♦ 0 while P(f) — ► Po > 0 
as £ — > 0 (see also Hosokawa and Yamamoto 1991). In other words, we test a 
question: what of the results from Sec. 2.2 which are obtained for e(x) are valid for 
the true dissipation £(x)? A comparison is executed in the following way. Using 
the DNS data base, each parameter is estimated for the true energy dissipation and 
one-dimensional surrogate by substituting in the same defining formula either e{x) 
or e(x), respectively. 
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FIGURE 11. DNS: pdd of the energy dissipation. In this and the following figures 
the solid symbols correspond to the true dissipation £(x), and the open symbols 
correspond to the one- dimensioned surrogate £(x). 




FIGURE 12. DNS: correlation functions of the energy dissipation, and the local 
values of scaling exponent /i r (fc). a , R\ = 94; A , 168. See caption of Fig. 11. 

The correlation functions R c (r) and the local values of scaling exponent /i r (r) are 
presented in Fig. 12. No clearly defined inertial range is seen, i.e., there is no region 
where ji r (r) % constant. This is expected due to relatively low Reynolds number in 
the present DNS. In spite of different pdd (Fig. 11), a qualitative behavior of iZ c (r) 
for £(x) and £(x) is quite similar. Hence one can expect that results and conclusions 
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FIGURE 13. DNS: pdd of the energy dissipation averaged over a segment r. 
(a) R\ = 94, r/r) = 10; (b) R\ = 168, r/r) = 15. See caption of Fig. 11. 



FIGURE 14. DNS: pdd of lne r . Solid lines correspond to the normal distribution, 
(a) R\ = 94, r/rj = 10; (b) R\ = 168, r/r} = 15. See caption of Fig. 11. 

of Sec. 2.2.1 are valid for the true energy dissipation field. 

Pdd of energy dissipation averaged over a segment r are presented in Fig. 13. 
Variation of i r (x) is much larger than that of e r (;r): the maxima of P{i r ) are shifted 
to the smaller amplitudes while the tails of P(i r ) at large amplitudes are higher 
than those of P(e r )- The log-normal model is tested in Fig. 14. The true energy 
dissipation field agrees quite well with the model while P(ln£ r ) reveals significant 
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r / r> r/r| 

Figure 15. DNS: variation of ln£ r , and the local values of /xjt(r) and (3(r). 
D , R\ = 94; a , 168. See caption of Fig. 11. 

deviations similar to those in Fig. 5. One can suggest that the deviations from the 
log-normal model in Fig. 5 would be smaller (or even disappear) if the true energy 
dissipation was measured at high R\. 

Eqs. (3) and (11) were tested, and the results are presented in Fig. 15 which is 
similar to Fig. 6 for experimental data. A qualitative behavior of <r, 2 n e for the one- 
dimensional surrogate in the DNS is completely similar to that in the high Reynolds 
number experiments. The DNS results for £ r do not agree with Eq. (3) at rjr) ss 
7—20 (this region is expected to model the inertial range) while they agree with 
Eq. (11). On the other hand, the true field e r (x) is reasonably described by Eq. (3) 
in the vicinity of rjr] as 10, and it does not agree with Eq. (11). Using DNS at 
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Figure 16. DNS: the non-centered and centered second-order moments of e r {x), 
and the local values of scaling exponents p £ (r) and 7e( r )- D > = A 1 

See caption of Fig. 11. 


R x = 68.1 and 151, Wang et al. (1994) recently found the values of Hk ~ 0.28 for 
the true field s r (x), and 0.34 for the surrogate e r (x). In the present DNS the values 
of y. k {r) for i r are also larger than those for e r (Fig. 15). However, our results differ 
conceptually from those by Wang et al. (1994). We state that the surrogate field 
does not obey Eq. (3) in the vicinity of r/r/ ~ 10-15, and no value of pr- may be 
estimated. Perhaps, Wang et al. (1994) would draw the same conclusion if they 
used the logarithmic derivatives instead of the log-linear plots of (r\ ntr us lnr. 

The DNS results for < e\ > (Fig. 16) are also similar to the experimental data 
(Figs. 8, 9). A functional behavior of the second-order moment < i 2 r > in the DNS 
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agrees with Eq. (13) at r/i) « 10-20 as that in experiments while <£ 2 > clearly 
agrees with Eq. (4). 

It follows from the above results that the experimental findings of Sec. 2.2.2 and 
2.2.3 may be valid only for (Du j D.v)~ , and their application to the true dissipation 
field s f r ) is quite questionable. Based on these preliminary results, we cannot draw 
any final conclusions. However, two suggestions may be stated. First, any exper- 
imental finding obtained for the surrogate i(x) should not be applied to the true 
dissipation field before similarity of this result to that for e(x) is proven. Second, 
DNS seems to be the most reliable (and perhaps a unique) tool for such a proof. 

3. Future plans 

The goal remains to understand what characteristics of the energy dissipation field 
are adequately represented by the one- dimensional surrogate. To further study this 
question, we are going to use the recent DNS data base with resolution 512 3 . 

The second question may be formulated as follows: what is the minimum set of 
terms in Eq. (2) which adequately represents any statistical characteristic of e(x )? 
An answer to this question will allow a proper design of the future high Reynolds 
number experiments. 
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